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^ ' Abstract 



We describe all the form factors for semileptonic and radiative B decays to the same 
light vector meson with just two parameters and the two form factors for semileptonic B 
decays to a light pseudoscalar meson with a further two or three parameters. This provides 
simple parametrizations of the form factors for the physical B^ —>■ p^l~ui, B^ Tr^l~i?i 
and B K*'-/ decays. The parametrizations are consistent with heavy quark symmetry, 
kinematic constraints and lattice results, which we use to determine the parameters. In 
addition, we test versions of the parametrizations consistent (or not) with light-cone sum 
rule scaling relations at = 0. 
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1 Motivation 



The aim of this note is to obtain a simple yet phenomenologically useful description of 
the form factors for semileptonic and rare radiative heavy-to-light meson decays for all g^, 
the squared four-momentum transfer to the leptons or photon. Lattice calculations provide 
values for the form factors over a limited region at high q^. We supplement these data with 
further constraints and model input to obtain our parametrization. 

The semileptonic decays have been measured by CLEO [|I| . Combining those results with 
our parametrization leads to values for \Vuh\- The decay 13 — >• K*-^ depends on \Vts\ when 
described in the Standard Model, but since it first occurs at one loop it is an excellent place 
to search for new physics 0, ^, ^ |^ : in either case it is important to know the relevant form 
factors. 

Kinematic constraints relate some of the form factors at = 0, and light cone sum rule 
(LCSR) calculations provide further information at low g^, particularly on the heavy meson 
mass dependence of the form factors as that mass becomes large. Specifically, LCSR predict 
that all the form factors scale like in the rriB ^ oo limit. At the other end of the range, 

near g^ax^ heavy quark symmetry (HQS) provides, at lowest order, extra relations between 
the form factors, but (unlike the case of heavy-to-heavy transitions) does not determine the 
overall normalisation. Extra assumptions are needed to cover the full range. For the 
decays with a light final state vector meson, we apply a simple model due to Stech , which 
uses just two free parameters. For a light final state pseudoscalar we use fits with two or 
three parameters 0, § . The normalisation in each case is determined using lattice data from 
the UKQCD collaboration g |, lig. 



2 Description of Model 

In this section, we transcribe the assumptions made by Stech in P] into the language of 
HQS. For q^ close to g^ax; ^^e standard leading order HQS analysis |]TT], [T2[ for a B meson 
decaying semileptonically or radiatively into a light pseudoscalar (vr) or vector (V) meson 
allows the corresponding matrix elements to be written as: 



{Ti\qi^h\B) = .JT^^Ti I \ei - ^V^7'| (1) 

and 

+ K + fe^r)^-^*]ri±^7^} , (2) 

where v is the four- velocity of the B meson, g is the light active quark {u or s), and F 
denotes either ^^{1 — 7^) or q^a^^{\ +7^)/2, with q^ = {ps —Pn,v)u- The 6^^'^ are functions 
of the invariant v ■ p^^.v ^^nd are independent of heavy-quark mass and spin, up to O (l/m^) 
corrections. 



7r,V 



In this language, Stech's model amounts to considering only the contribution from 62 
and setting all other 9's to zero. In this sense, the model is a truncation of the HQS results. 
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As we shall see, its usefulness and justification rest on the facts that it satisfies all known 
constraints and fits the lattice data (available at large q^) well. Evaluating the traces in 
Eqs. (|1],0) and solving the resulting equations for the form factors in a standard notation, 
we find 

Ao{q') 
V{q') 
Ai{q') 

A2{q') 
2T,{q') 
2iT2{q^) 

which are equivalent to the relations given in Eq.(ll) of Ref. with, in addition, relations 
involving Ti and T2. 

As we will see below, these relations can be used for a simple and successful description 
of lattice results for the form factors. We recall that, by construction, this model respects 
HQS scaling relations and satisfies the two kinematic constraints: 

Fi(0)=Fo(0), Ti{0)=tT2{0) (4) 

For the parametrization to be complete we must specify the functions ^2'^. This will be 
discussed in Section |. 

3 Lattice Details 

The details of the lattice simulation can be found in Ref. ^j. The chiral extrapolations for 
— > p^l~iii and B K*'y are described in Ref. while those for B^ are 
explained in Ref. 0. The only way in which our determination of the lattice form factors 
differs here is in the heavy-quark mass extrapolation. We improve this extrapolation for the 
form factors Ai, V, Ti and T2 by imposing heavy quark symmetry in the infinite heavy-quark 
mass limit, as we now explain. 

All form factors are calculated for four values of the heavy quark mass around the charm 
mass and for a variety of values of g^. In previous work 0, the form factors were extrap- 
olated at fixed four- velocity recoil, uj = v ■ {p-K,v/'m-Ky)i near the zero recoil point u; = 1, by 
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fitting to the following heavy-quark scaling relations: 

M ^ M2 



/eM-//^ = 7/(i + ^ + ^ + ---) (5) 



where / = Fi, Fq, Aq, V, Ai, A2, Ti, T2 and Nf = — 1, 1, — 1, — 1, 1, — 1, — 1, 1 respectively. 
Here, M is the mass of the heavy-light meson, 7/, 6f and e/ are the fit parameters^, and G 
comes from leading logarithmic matching and is chosen to be 1 at the B mass ||13|| , 

with /3o = 11 in the quenched approximation and Aqcd = 200 MeV. 

While correct, this procedure neglects the fact that in the M — > cxd limit, HQS predicts 

A,iq')=2tT,{q'), V (q') = 2T,iq') (7) 

for q^ not too far from g^ax; i-^- close to zero recoil. In the present paper, we enforce this 
prediction of HQS by performing a combined fit, at fixed u, of the pairs of form factors 
{Ai,T2) and {V,Ti) to the parametrizations of Eq. with the following constraints on the 
fit parameters: 

7Ai = 2i7T2, 7v = 27ti • (8) 

This not only guarantees that the extrapolated form factors are consistent with HQS in the 
infinite mass limit, but also reduces statistical errors, because the number of parameters is 
decreased. 

As an example, we show in Fig. ^the combined fit of (V,2Ti) to the constrained heavy 
quark scaling relation just described, for the case of a final state p. Both the linear and 
quadratic fits in 1/M are excellent, confirming our earlier finding, in Ref. |^, that the form 
factors Ai, V, Ti and T2 satisfy the infinite mass relations of Eq. (|^) very well, even when 
extrapolated independently. Because linear and quadratic extrapolations always give results 
which agree within error, because the quadratic term is always consistent with and because 
the quadratic fits may emphasize small discretization errors, we will use the linear results in 
the following. 

In the results reported below, all errors are statistical central 68% bounds. For decays to 
a light final state vector we use 250 bootstrap samples from our lattice data. For decays to 
a light pseudoscalar we add errors to account for the chiral extrapolation 0] and propagate 
them assuming standard Gaussian statistics. 



4 Form Factors 

One could imagine performing a combined fit to n^l^ui, B^ p^l^vi and B K*'~f 

form factors assuming ^J(5^)=^^(5'^) in Eq. (^). However, there is no reason to have spin 
symmetry relating the final pseudoscalar and vector states. Furthermore, our lattice results 
clearly indicate that the q^ dependences of Fi{q^) and Aoi^q^) are different, in contradiction 
to the assumption Ol^i^q^) = O'^^q'^). Thus we consider decays to pseudoscalar and vector 
states separately. 

^e/ is set to for linear extrapolations. 
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Figure 1: Constrained l/M extrapolation of the pair of form factors {V,2Ti), as de- 
scribed in the text. The sohd curves correspond to a hnear extrapolation, the dashed 
curves to a quadratic fit. The four lower circle points are the lattice results for 
2Ti^JmB/M and the four upper ones are those for V^JmB/M, in the case where a B 
meson at rest decays to produce a p meson with one (lattice) unit of three momentum. 
The points at nip/M = 0.146 correspond to the linearly (diamonds) and quadratically 
(squares) extrapolated results for the two form factors at M = tub- 
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0.04 U.ZO_o.03 U.i)U_o.04 


0.35± 


0.06 
0.05 


g2=0 0.16± 


0.02 
0.01 


?max 


0.46± 


0.02 n 88+0-05 1 on+0.09 
0.01 U.08__o.03 -L-oU„o.05 


2.071 


0.11 
0.06 


gLx 0.90± 


0.05 n 9c;+ 0.01 

0.04 U.ZO_o.01 



Table 1: Results of fits to tfie lattice predictions for Aq, Ai, V, Ti and T2 assuming a 
pole form for Ai, as given in Eq. (|^), in the parametrization of Eq. (^). On the left are 
results for a p meson final state (Situation A) and on the right, results for a K* meson 
final state (Situation B). We also provide the values of the form factors for p'^l~vi 
and B — > K*'^ decays at g^=0 and ^max; ^ determined by the fit. 



4.1 ^ p^l vi and B K*^ decays 

So as not to assume flavour S't/(3) symmetry in our combined description of p^l^vi 
and B K*'y decays we have used the freedom to adjust quark masses in lattice calculations 
and considered two situations: 

A The mass of the light active quark q in Eq. (^ is set to the u quark mass (i.e. V = p 
in Eqs. (@,|D). In this case, the combined fit to the form factors of the operators 
g7'^( 1—7^)6 and gg^o"^*^ (1+7^)6/2 constrains, through Eq. (^), the form factors relevant 
for the decay B^ p^l~ui. 

B The mass of the light active quark q in Eq. (^ is set to the s quark mass (i.e. V = K* 
in Eqs. (H,^). Now the combined fit constrains, through Eq. (|^), the form factors 
relevant for the decay B — > K*^. 

For the parametrization of Eq. (|^) to be complete we must specify the function 62 or, 
equivalently, one of the form factors. With pole dominance ideas in mind, we shall consider 
a pole form for Ai(g^) which is compatible with previous lattice analyses ||^ and LCSR 
results fT^ : 

M<l') = , (9) 



where the two free parameters are ^i(O) and Mi, a mass on the order of niB- With these 
two parameters alone, we can describe the six form factors needed for decays to a vector 
final stateQ. This form further guarantees that all form factors scale like vn^^'^ at g^ = 
as predicted by LCSR in Ref. |T^. The results for 5° p^l~i^i and B K*'y decays are 
presented in Tab. |I], and Figs. |^ and ^ As the indicates, this simple parametrization works 
surprisingly well. The SU{3) flavour dependence in our results, from comparing situations 
A and B, is mild. Corresponding form factor values at g^ = differ by less than 10%. 



^In practice, we fit lattice data for the five form factors V, Aq, Ai, Ti and T2 to fix our free parameters. 
Our present lattice results do not allow a reliable extraction of A2 . 
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Figure 2: Fit to the lattice predictions for Aq, Ai, V, Ti and T2 for a p meson final state 
(Situation A) using the parametrization of Eq. (^) and assuming a pole form for Ai as 
given in Eq. @. The dashed vertical line indicates Qmax- 
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Figure 3: Fit to the lattice predictions for Aq, A\^ V, Ti and T2 for a K* meson final 
state (Situation B) using the parametrization of Eq. (|3|) and assuming a pole form for Ai 
as given in Eq. (^). The dashed vertical line indicates ^^ax- 
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5° ^ p+l-ui B K*-f 

Ao{0) = OAOt^ol Ao{0) = QA2±rol 

Mo = 5.10± Hi GeV Mq = 5.05± GeV 
xVdof = 25/20 xVdof = 27/20 



A2 Ao V Ti T2 



q'=0 0.35±g:gi 0.34±0:gi 0.40±g:g| 0.45±g:g| q'=0 0.21±g:01 

gLx 0.46m? 0.89±m 1.81±U', 2.08±Ul qL. 0.91±^t 0.251^1 



Table 2: Results of fits to the lattice predictions for Aq, Ai, V, Ti and T2 assuming 
a pole form for Aq, as given in Eq. (0), in the parametrization of Eq. (|^. On the left 
are results for a p meson final state (Situation A) and on the right, results for a K* 
meson final state (Situation B). We also provide the values of the form factors relevant 
for p'^l^h'i and B K*^ decays at q^=Q and g^ax determined by the fit. 



We have considered functional dependences for Ai different from that in Eq. (^. Constant 
behaviour has aheady been ruled out by lattice results Ref. Dipole (and in general higher 
powers: tripole, . . . ) and pole fits are hardly distinguishable in the physical range < < 
g^a^. In a dipole fit, the mass parameter Mdipoie is roughly given by Mdipoie ~ v^^poic and 
hence pole and dipole fits agree in the range of values of explored in Figs. ^ and |^. We 
have also studied "modified-constant" behaviour for Ai. 

A{q') = {l-qVml) "^f) , , (10) 

1 — q'^/m^ 

where rua and rri/s are masses on the order of ms. Such a behaviour was introduced in Ref. 
in the context of B ^ irlu decays. Within the framework of Eq. (^), this can be implemented 
by parametrizing Ao[q^) as a pole without increasing the number of free parameters. Thus, 
we have performed a second set of fits where we complete the parametrization of Eq. (|^) by 
taking the following functional form for Aq: 

Mi) - (11) 

where the two free parameters are Ao(0) and Mq, a mass of the order of m^. The cor- 
responding results for B^ p^l'J^i and B K*'y decays are presented in Tab. ^. This 
parametrization yields different behaviour for the form factors at small q^: the values at 
q^ = differ by two to three standard deviations. 

Though we cannot discriminate against this parametrization on the basis of alone, 
there are physical arguments for preferring the Ai pole fits. Fixing pole behaviour for Aq in 
the context of Stech's parametrization forces all form factors to diverge at the same value 
of even though different form factors receive contributions from resonances with different 
quantum numbers and masses. Furthermore, we find that this divergence occurs at around 
5.1 GeV, below the first physical pole, the B. Assuming pole behaviour for Ai, on the other 
hand, allows the form factors which receive contributions from 1+ resonances {Ai, A2 and 
T2) to diverge at larger values of q^ than the more singular form factors V , Ti and Aq. 
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Our fit confirms such behaviour as we find Mi = 7.0± GeV, in reasonable agreement with 



lattice [|T5l and quark model estimates [jT6|, |T^, [T^ for the l"*", bq resonance. The more singular 



form factors then diverge at = + my which is close to the physical 0~ (Aq) and 1~ (V^ 
and Ti) poles. Finally, LCSR scaling relations at = combined with HQS requirements 
at g^a^x ^i^l^ out pole behaviour for Aq but not for Ai. Thus, for physical applications, we 
consider only the parametrization given in Eqs. (|^J^). 

4.2 Ti^l'vi decays 

Stech's parametrization predicts that -fo(?max) vanishes in the chiral limit, in contradiction 
with our results and made unlikely by unitarity bounds |^. Furthermore, the B* which 
contributes a pole to Fi induces the same singularity in Fq within Stech's model. Because 
the B* pole is very close to g^ax' provokes a much more pronounced g^ dependence for Fq 
than observed in the lattice results or induced by the nearest resonance in this O"*" channel 
whose mass is significantly larger than itlb*- 

Restricting to polar-type g^-dependences, consistent with the kinematical constraint, 
Fi(0) = -Fo(O), heavy quark symmetry and unitarity bounds [0, H, we consider two pos- 
sible functional forms: 



pole/dipole 



(1 - qymf) l-L - 9 /"^oJ 



"modified-constant" / pole 

= 7T^^^ ^ Fo{q') = i^(0) (13) 



(1 — q'^/mfj (1 — g^/mo 



These two dependences have been studied previously in Ref. with light quark masses 
slightly larger than that of the strange quark and in Ref. for massless u and d quarks. The 
fit results are summarized here in Tab. |^. We note that only the pole-dipole parametrization 
is consistent with LCSR scaling relations at g^ = together with HQS requirements at g^ax- 
Both fits are consistent with the dispersive bounds of Ref. M. 



5 Phenomenological Consequences 

Using the preferred pole form for Ai in the decay B'^ p'^l~ui, we find the differential decay 
rate spectra in g^ and the lepton energy E.^ These are shown in Fig. ^ Likewise, using the 
pole/dipole fits for 5° — vr we find the differential spectra shown in Fig. |^. Integrating 
to find the total decay rates we have, for massless leptons in the final state: 

T{B'^p+l-ui)/\V^,\^ = 10.9±?:ix 10-12 GeV =16.5±|ips-i (14) 
r(S° ^ 7r+/-z70/|Kfer = 5.6lg:ix 10-12 GeV = 8.5±ijips-i (15) 

^The errors quoted here differ from those in which were determined by an incorrect procedure. The 
central fit values agree. This difference does not affect the dispersive bounds of Ref. 0, nor the comparison 
of these bounds with the various parametrizations. 

*See Ref. Il^ for the decay rate kinematics, including the case of non-zero lepton masses. 
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fit type 



F(0) mi(GeV) mo ( GeV) ma ( GeV) xV^of 




5 10 15 20 



g2(GeV') 




E{ GeV) 



"modified-const" 0.43 ± 0.06 m^* 5.5 5.89 ± 0.52 0.5/4 

pole/dipole 0.27 ±0.11 5.79 ± 0.58 6.1 ± 1.5 - 0.1/3 

Table 3: Fits of the lattice results for Fi{q^) and FQ{q^) to the parametrizations described 
in the text (taken from Ref. Q). In order to perform the "modified-constant/pole" fit 
with a limited range of results, mi and rriQ were fixed to the value of the corresponding 
nearest resonances. (This latter fit is named "fixed-pole" in Ref. 0.) 



Figure 4: Differential decay spectra for p'^l vi for massless leptons: (a) dV/dq^ 

in units of 10-i2|14^|2 GeV-\ (b) dV/dE in units of W-^'^\Vub\'^ . The form factors are 
taken from the fit to the parametrization of Eq. (^, assuming a pole form for Ai as given 
in Eq. (^. The dashed lines show the envelope of the 68% bootstrap errors computed 
separately for each value of q'^ or E respectively. 



For decays with a tau lepton in the final state we find: 

r(5° ^ p+r-z7,)/|K6r = 5.8±[j|x 10-12 = 8.81 Jips-^ (16) 
r(S° ^ n+r-iyr)/\Vub\^ = 3.8±li x lO^^' GeV = 5.8i ^ips^^ (17) 



In Tab. ^ we compare our results with recent LCSR |T^, ^ ^ |2^ calculations^. 

The agreement of the form factor values at = is very good. In the spectrum dr{B^ —>■ 
p^l^ui) / dq^ we note that our result is larger near = 0, as is our value for the ratio Tl/Tt- 
The differential decay rate at = is proportional to 

which is the longitudinal contribution, so that small differences in Ai{0) and ^2(0) can lead 
to significant differences in the differential rate and in the ratio of longitudinal to transverse 



^We quote LCSR results at leading order in perturbative QCD. The O (as) corrections to the leading 
twist term only are now known for Fi |25| . 
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g2(GeV2) E{GeV) 

Figure 5: Differential decay spectra for tt^I^ui for massless leptons: (a) dV /dq^ 

in units of lQ-^'^\Vub? G<^~^ , (b) dV/dE in units of lQ-^'^\Vub? ■ The form factors are 
taken from the fit with a dipole form for Fi as given in Eq. ([l^). The dashed Unes show 
the envelope of the 68% errors computed separately for each value of or E respectively. 



rates. Since the lattice results are available for large g^, the dV /dq^ plots in Figs. | and ^ 
also show errors growing larger towards = 0. This is complementary behaviour to LCSR 
results, which are more reliable at low q^. 

All other rates agree with the LCSR results within errors. Agreement is also excellent with 
the experimental results of CLEO For the ratio r(p)/r(7r), which is independent of \Vub\ 
and can therefore be compared directly to the value in Tab. |, CLEO quotes 1.4±o;4±0.3±0.4, 
where the errors are statistical, systematic and the estimated model dependence. Moreover, 
branching ratios, obtained from the rates in Tab. § with the values of \Vub\ and of the B 
lifetimes given in Ref. |1| by CLEO, compare very favorably with the values in that same 
publication. 

Finally, the additional constraints provided by our approach may shed some light on 
the twofold ambiguity in the lattice value of Ti(0) (see for example |2^) by favouring 
the smaller of the two determinations. With our prediction for Ti(0), we can determine 
the hadronization ratio Rk* = r(i? — >■ K*'j)/r{b — > 57), which is given, up to 0{l/ml) 
corrections, by 

fl,. = 4(!^iiiyfi-!%)'|r.(0)r (19) 

We find Rk* = (1613)%, to be compared with the experimental value (18 ± 7)% p8| . 

6 Conclusion 

Inspired by the work of Stech [Q, we have designed a simple parametrization for the form 
factors which describe B transitions, where is a light vector meson. This parametriza- 
tion is consistent with heavy quark symmetry and kinematical constraints but requires an 
ansatz for the g^-dependence of one of the form factors. The parameters of the ansatz are 
determined by fitting to lattice results around g^ax- We have explored several ansatze and 
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Fi(0) 



Am 



^2(0) 



Ao(0) V{0) 



ri(o) 



This work 
LCSR 
LCSR 

LCSR im, m 

LCSR g 



0.27 ±0.11 



0.24-0.29 



n 97+0.05 

0.27 ±0.05 
0.24 ±0.04 



n 9(3+0.05 
0.03 

0.28 ±0.05 



0.30±g:gt 



n oc+O.Oe 
U-JO- 0.05 

0.35 ±0.07 
0.28 ±0.06 



0.16: 



-0.02 
0.01 



0.16 ±0.03 



0.15 ±0.03 



T 



This work 
LCSR 



8.5: 



3.3 
1.4 



3.5 
2.3 



m 



16.5^ 
13.5 ±4.0 



1.7 ±0.5 



0.80±^i 
0.52 ±0.0^ 



LCSR 231 8.7 



Table 4: Form factor values at = and decay rates and ratios for 6 — > u transitions 
from this calculation and from light cone sum rules (LCSR). Decay rates are given in 
units of |T4f,p ps""*^. 



favour pole behaviour for Ai. though not singled out by alone, it satisfies LCSR scaling 
relations at = and allows form factors that receive contributions from resonances with 
different quantum numbers and masses to diverge at different g^. As a result we describe 
semileptonic and radiative B decays to the same light vector meson with only two param- 
eters. We use the freedom provided by lattice calculations to consider two independent 
situations. In the first, the mass of the final state vector meson is set to the physical p mass. 
We then apply our model and fit both the semileptonic and radiative form factors simul- 
taneously, keeping the former to describe the physical B^ p^l~vi decay. In the second 
situation, the final state meson is chosen to be a K* . Again we fit all form factors at once, 
this time keeping the radiative form factor results to describe the physical B — > K*'^ decay. 
In this way, we do not assume S'f/(3) flavour symmetry. 

For semileptonic B transitions to a light pseudoscalar meson, we perform a separate 
analysis, since we do not assume spin symmetry for the final state meson. Our preferred 
description, consistent with heavy quark symmetry, kinematic constraints and LCSR scaling 
relations, is a pole/dipole fit with three parameters. We apply this to B^ ti^I~i)i decays. 

Even though our approach requires some assumptions, we have tried to minimize their 
number and make them consistent with all known theoretical constraints and lattice results. 
The resulting parametrizations should be extremely useful for phenomenological applications 
because they provide a simple and effective description of the behaviour of the various 
form factors. 

We should also like to stress that we have introduced an improved procedure for extrap- 
olating form factors in heavy-quark mass from around the charm, where they are calculated 
on the lattice, to the 6-quark mass. This procedure reduces statistical errors by making full 
use of the constraints of heavy quark symmetry. 
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